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Motivation 

Many progressive failure analysis (PFA) tools have been 
developed and used; however, predicting the failure 
load of composite structures is still a challenging problem 

Since damage initiation in composite laminates often involves 
out-of-plane stresses such as the familiar cases of edge 
delamination and low velocity impact damage, 3D 
progressive failure analysis tools are highly desired 

Recent studies 1 show that lamina property degradations 
after damage initiation can be modeled by the crack band 
theory 2 that can also alleviate mesh dependency issues 

Crack band theory provides the key ingredient for establishing 
the continuum damage mechanics (CDM) based progressive 
failure tool 


1 .Pineda and Waas, Int J Fract (201 3) 
2. Bazant and Oh, Mater Struct (1983) 


Objective 


Develop a CDM based 3D PFAtool that can be linked 
to commercial finite element analysis (FEA) software 
for laminated composite structures 


Approach 

Develop CDM model based 3D constitutive equations 
for damaged lamina, utilizing only standard unidirectional 
material properties 

Implement the CDM model as a user defined material 
model (VUMAT) to link with Abaqus/Explicit code for 
performing PFA of laminated composite structures 

Validate the PFA tool with available experimental data 


CDM Based PFA 


Strain based failure criteria for damage initiation 

Damage evolution modeling 12 

- Crack-band strain in damaged element 

- Crack band theory for material property degradations 

Matzenmiller type 3D stress-strain relations 1 2 3 
for damaged lamina 


1 . Pineda and Waas, Int J Fract (2013) 

2. Bazant and Oh, Mater Struct (1 983) 

3. Matzenmiller, Lubliner, and Taylor, Mech Mater(1995) 


Strain Based Failure Criteria for Damage 

Initiation 

• Hashin-Rotem Failure Criteria 
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Crack-Band Strains in Damaged Element 12 


Smearing the crack opening displacements to 
increase the element strains, namely, the crack- 
band strains 
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Damaged element with cohesive opening 

l e = length scale for ply in - plane strains 
V = length scale for ply transverse strains 


Crack-band strains combined with the traction- 
separation laws for determining material property 
degradations 


1 . Pineda and Waas, Int J Fract (2013) 

2. Bazant and Oh, Mater Struct (1983) 


Crack Band Theory for Fiber Direction 

Modulus Degradation 


Linear elastic Damage evolution 
modulus / , 



Crack-band strain 
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Traction-separation law 



Material Moduli Degradations 
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Pineda and Waas, Int J Fract (201 3) 


Damage Indices (Internal State Variables) 


Definition of damage index (Internal state Variable) 

d =\-E IE 
u i 1 Mi ' Mio 


^2 1 ^22 f ^220 

d 3 —1 — E 33 / £ , 330 


^12 — 1 _ ^12 ^ ^120 
^13 = 1 _ ^13 ^ ^130 
= 1 _ ^23 ^ ^230 


Zs n , E n 2 , £33, G p , G 13 , a/ir/ G 23 are degraded material moduli 

E n 0 , £' 22() , Me) > Gpo? G 130 , aad G 230 are linear elastic material moduli 



• Note that damage index = 0 is at undamaged state (fully elastic) 

damage index = 1 is at fully failed state 

• Continuum damage mechanics model using these damage 
indices to establish the material stiffness degradations 


Stress-Strain Relations for Damaged Lamina 


o = C(d) e where C(d) 



C 

C 


12 


22 


c 


13 


c 

c 


23 


33 


Symmetry 


C n = (1- d^Efi - (1- d 2 )(l- d 3 )v 32 v 23 ] / A 
C 2 2 = (1 — d 2 )^2 [1 — (1 — d l )(1 - d 3 )v 13 u 31 ] / A 


C 44 = 2d 4 G l2 / A 
C 66 =2d 6 G 23 /A 


0 0 0 

0 0 0 

0 0 0 

c 44 0 0 

C 55 0 

C 66_ 

C 55 =2d 5 G 3l /A 


C 33 = (1 -d 3 )E 3 [1 - (1 - d { )(1 - d 2 ) v 12 o 21 ] / A 
C 12 = (1 — )(1 — <^ 2 )£’ 1 [(1 — <^ 3 ) V 31 \3 23 + u 21 ]/ A 

C 13 = (\ — d l )(\ — d 3 )E x \(X — d 2 )\ 2 {0 32 +o 31 ]/A 
C 2 3 = (l- d 2 )(l- d 3 )E 2 [(l- d x )\ n u 3l +o 32 ]/A 

A = 1 — (1 — d 2 )(l-d 3 )v 32 v 23 -(1- d { )(\- d 3 )v 3l v l3 -(1- d x )(l- J 2 )v 21 v 12 
- 2(1 - d x ){ 1 - J 2 )(l - d 3 )v 13 v 21 v 32 


d 4 ,d 5 , and d 6 can be expressed as the 
following phenomenological equations 

d 4 = (l-J 1 )(l-oui 2 )(l-pds 12 ) 
d 5 = (1- J 1 )(l-ad 3 )(l-pds 13 ) 
d 6 = (l-J 1 )(l-aJ 2 )(l-ad 3 )(l-pds 23 ) 


Matzenmiller, Lubliner, and Taylor, Mech Mater(1995) 


VUMAT for Abaqus/Explicit Analysis 


Implementing the CDM based 3D PFA equations 
developed as a user defined material model (VUMAT) 
in Abaqus/Explicit Code 

Implemented for the reduced integration 3D element, 
C3D8R 

Applicable for geometrically nonlinear analysis 

Using mass scaling to shorten the analysis time while 
controlling the proper time incrementation for 
obtaining stable accurate solutions 


Validation of Analysis Results 


AFRL open hole specimens test data 

- From publications 1-7 related to AFRL program entitled 
“Damage Tolerance Design Principles (DTDP)” 


AIAA SciTech 2015 (56 th SDM conference) 

1 . Clay and Holzwarth 

2. Engelstad, Action, Clay, Holzwarth, Robbins, and Dalgarno 

3. Zhang, Patel, and Waas 

4. Joseph, Waas, Ji, Pineda, Liguore and Wanthal 

5. larve, Hoos, Braginsky, Zhou, and Mollenhauer 

6. Fang, Cui, and Luya 

7. Abdi, Codines, DorMohammadi, and Minnetyan 


IM7/977-3 OHT and OHC Specimen 


138 mm 


0-deg ply 

38.1 mm ^ ) 

6.35 mm 

y 

Clamped 

Edge 



Three different layups 

1 . [0/45/90/-45] 2S 

2. [60/0/-60] 3S 

3. [30/60/90/-60/-30] 2S 


Clay and Holzwarth, AlAASciTech Conference 2015 


3D Finite Element Mesh of Open Hole 

Specimen 


u 



• One element per ply thickness 

• Only top half of the laminate modeled (symmetry boundary conditions employed) 


IM7/977-3 Material Properties 12 


Material 

Property 


E u (GPa) 

164/137.4 

E 01 ,E 33 (GPa ) 

8.977 

G l2 , G l3 (GPci) 

5.02 

G 13 (GPa) 

3.00 

Vl2>V 13 

0.32 

V 23 

0.496 

G IC , G IIC (N / mm) 

0.256, 1.156 

Y T ,Y c (MPa) 

100.0, 247.0 

S(MPa) 

80.0 

X T ,X c (GPa) 

2.9, 1.68 

G FI , G fc ( N / mm) 

81.534, 24.533 

Gmt ’ G mc ! mm) 

0.256, 1.156 


1. g ft ,g fc : Fiber tension and 
compression toughness parameters 

2. G mt ,G mc : Matrix tension and 
compression toughness parameters 


1 AIAA SciTech 201 5 Papers: 

a) Clay and Holzwarth 

b) larve, Hoos, Braginsky, and Zhou 

c) Zhang, Patel, and Waas 

2.Camanho, Mami, Davila, Compos Sci 
Technol, (2007) 


Stress-Strain Curves for Various Beta Values 

([0/45/90/-45] 2S OHT) 



0 0.25 0.5 0.75 1 1.25 1.5 

Strain (%) 


Same initial stiffness for all 
cases 

As p increases, predicted 
failure stress decreases 
All predicted failure stresses 
exceed 90% of the test 
maximum stress 


Using the distance of extensometer knife edges, located half an inch 

(along the length direction) below and above the center of the hole for computing 

strain 


Test and Analysis Correlation of Maximum 

Applied Stress 

Performed many parametric studies with various values of 
a and p and selected a = p =0.5 for shear stiffness degradations 
for all OHT and OHC cases analyzed 


Lay-Ups 

OHT 

OHC 

[0/45/90/-45] 2S 

-8.82% 

-6.94% 

[60/0/-60] 3S 

-2.52% 

-18.3% 

[30/60/90/-60/-30] 2S 

6.03% 

8.32% 


Reasonable good correlations obtained; however, further development 
or adoption of more accurate PFA tool needed 


Influence of Beta on Damage Evolution in 

0-Plies ( [0/45/90/-45] 2S ) 



Assume 
a = p 


Load and 0-ply 
direction 


At 82% maximum stress 


Outer 0-Ply 


d, index 



+ 1.000e+00 
+7.500e-01 
+ 5.000e-01 
+2.500e-01 
+0.000e+00 


Inner 0-Ply 

P = 0.2 "^p = 0.5 

As p increases, damage progression increases for same load level 





Concluding Remarks 

Developed and implemented a CDM based explicit 3D PFAtool as 
a VUMAT subroutine to work with 3D elements in a commercial 
FEM code 

Key features implemented in the VUMAT including: 

- Simple failure criteria used for failure initiation 

- Lamina material properties obtained from standard test procedures 

- Technique for alleviating mesh dependency implemented 

- Damage due to 3D stress fields considered 

- Matzenmiller type stress-strain relations for damaged lamina 

used in which the stiffness degradations expressed 
as functions of the fiber, matrix, and shear damage indices 

Preliminary predicted failure stresses having reasonably good 
agreement with test data of OHT and OHC specimens 
with various layups 


